Observation of the Decay Kl fJ.~^fJ. 77. 
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We have observed the decay Kl — > /x"^^~77 at the KTeV experiment at Fermilab. This decay 
presents a formidable background to the search for new physics in Kl —* tv''^ fi'^ fi~ . The 1997 data 
yielded a sample of 4 signal events, with an expected background of 0.155 ± 0.081 events. The 
branching ratio is B{Kl fi'^n'yy) = (10.4l5;g (stat) ± 0.7 (sys)) x 10"*^ with m^-y > 1 MeV/c^, 
consistent with a QED calculation which predicts (9.1 ± 0.8) x 10~^. 



In this paper we present the first measurement of the branching ratio for Kj^ — > n'^ ij,~"f"f. This decay is expected 
to proceed mainly via the Dalitz decay — > fj.'^ fi^^ with an internal brcmsstrahlung photon. This decay is one of a 
family of radiative decays {Kl /i+/i^7, Kl /i+/i^77, Kl e+e^7, Kl — ^ 6+6^77) which are under study at 
KTeV and elsewhere 0J^. The decay Kl fi'^fi^^^ presents a formidable background to the search for direct GP 
violation and new physics in Kl — > n^fi^fj," decays 

The measurement presented here was performed as part of the KTeV experiment, which has been described else- 
where 0]. The experiment used two nearly parallel Kl beams created by 800 GeV protons incident on a BeO target. 
The decays used in our studies were collected in a region approximately 65 meters long, situated 94 meters from the 
production target. The fiducial volume was surrounded by a photon veto system used to reject events in which photons 
missed the calorimeter. The charged particles were detected by four drift chambers, each consisting of one horizontal 
and one vertical pair of planes, with typical resolution of 70 fim per plane pair. Two drift chambers were situated on 
either side of an analysis magnet which imparted 205 MeV/c of transverse momentum to charged particles. The drift 
chambers were followed by a trigger hodoscope bank, and a 3100 element pure Csl calorimeter with electromagnetic 
energy resolution of a{E)/E — 0.45% © 2.0%/ EiGcV). The calorimeter was followed by a muon filter composed 
of a 10 cm thick lead wall and three steel walls totalling 511 cm. Two planes of scintillators situated after the third 
steel wall served to identify muons. The planes had 15 cm segmentation, one horizontal, the other vertical. 

The trigger for the signal events required hits in the upstream drift chambers consistent with two tracks, as well as 
two hits in the trigger hodoscopes. The calorimeter was required to have at least one cluster with over 1 GeV in energy, 
within a narrow (20 ns) time gate. The muon counters were required to have at least two hits. In addition, preliminary 
online identification of these decays required reconstruction of two track candidates originating from a loosely-defined 
vertex, and each of those track candidates was required to point to a cluster in the calorimeter with energy less than 
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5 GeV. A separate trigger was used to collect Kl 7r+7r~7r'' decays which were used for normalization. This trigger 
was similar to the signal trigger but had no requirements on hits in the muon hodoscopes or clusters in the calorimeter. 
The preliminary online identification was performed on the normalization sample as well, but no energy requirements 
were made on clusters pointed to by the tracks. The normalization mode trigger was prescaled by a factor of 500:1. 

The main background to Kl — > fJ,'^ fi~"f^ was the Dalitz decay — > /i^/x~7 with an additional cluster in the 
calorimeter coincident with but not due to the decay. Such an "accidental" cluster could appear as a photon. 
Additional backgrounds were Tr+vr^vr*' decays with the charged pions misidcntified as muons due to pion decay 

or pion punchthrough the filter steel, and — > n^fj.^i' decays {K^^) with both charged pion misidentification and 
accidental cluster contributions. Other contributions, such as Kl tt+tt^ decays and Kl 7r+7r^7 decays, were 
negligible. 

Offline analysis of the signal required the full reconstruction of exactly two tracks. The vertex reconstructed from 
the two tracks was required to fall between 100 meters and 158 meters from the target. In order to reduce backgrounds 
due to pion decay in flight, we required that the track segments upstream and downstream of the analysis magnet 
matched to within 1 mm at the magnet bend plane. Further, we required the calculated from the reconstructed 
two-track vertex be less than 10 for 1 degree of freedom. Tracks were required to have momenta equal to or greater 
than 10 GeV/c to put them above threshold for passing through the filter steel but below 100 GeV/c to ensure well 
measured track momenta. Since muons typically deposit ~ 400 MeV in the calorimeter, we required the energy 
deposited by each track be 1 GeV or less. In addition, we required two non-adjacent hits in both the vertically and 
horizontally segmented muon counters. 

Figure |l| shows the expected distribution of cluster energy due to photons from Kl — ^ events and those 

from accidental sources. Accidental clusters in the calorimeter were typically of low energy. Events were required 
to have two calorimeter clusters consistent with photons with no tracks pointing to them. One of these clusters was 
required to have greater than 10 GeV of energy, thus reducing backgrounds due to accidental clusters. 



In order to reject backgrounds from decays that contained a tt , the invariant mass of the two photons, m 
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was required to be less than 130 MeV/c^. Approximately 8% of the Kl 7r+7r~7r° decays in which the charged 
pions decay to muons survived the m-yj cut because the mismeasurement of the charged vertex smeared the m-yj 
distribution. In order to remove these events, we constructed a variable (Rjj''^) defined as 

where rriK is the kaon mass, mTTTr is the invariant mass of the two tracks assuming they are due to charged pions, 
p'±-KTr square of the transverse momentum of the two pions with respect to a line connecting the target to the 

two-track vertex, and uIt^o is the mass of the tt^. This quantity is proportional to the square of the longitudinal 
momentum of the ir^ in a frame along the Kl flight direction where the tt+tt" pair has no longitudinal momentum. 
Figure |2| shows the expected distribution of KJ^'^ for the signal (using an 0{a) QED matrix element), and the 

Kl — > 7r+7r~7r*' background. By requiring R||'^ to be -0.06 or less, 92.7% of the remaining Kl Tr+7r^7r° background 
was eliminated. 

The invariant mass of the two tracks assuming muons, m^^, provided a way to reduce backgrounds due to K^^ 
decays. Figure ^ shows the expected distribution of for the signal and background. We required to^^ to be less 
than 340 MeV/c^. This cut ehminated 92.9% of the iC^g events. 

The cosine of the angle between the two photons in the kaon rest frame, cos6^~f, was also used to reject K^j^^ decays. 
The distribution of cosOj^ for the signal peaks at -1 corresponding to anti-coUinear emission of the two photons. The 
Kfj_3 background, which has two accidental clusters identified as photons, displays no such correlation. Figure ^ shows 
the expected cos6'^^ distribution for signal and K^^ background. We required cosO^^ to be -0.3 or less. This cut 
rejected 85.3% of the remaining Kfj_3 events. 

We also required the transverse shower shape for the photon clusters to be consistent with that expected from an 
electromagnetic process. The of the spatial distribution of energy deposited in the calorimeter was used to identify 
clusters as photons. This cut reduced the remaining backgrounds due to accidental energy by a factor of 4.5 while 
retaining 98.8% of the signal events. 

In order to estimate the amount of background in the signal region, we simulated all the leading sources of back- 
ground. Our simulation incorporated both charged pion decay in flight and punch-through the fllter steel. The 
punch-through rate was a function of tt^ momentum, determined by a Kl — * T^ei' control sample. The effect of 
accidental activity was simulated by overlaying Monte Carlo events with data from a random trigger that had a rate 
proportional to the beam intensity. The estimated background level is detailed in Table ||. A total of 0.155 ± 0.081 
background events are expected within the signal region. This region is defined by the invariant mass of the /i"*'/i~77 
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(m), and square of the transverse momentum of the four particles with respect to a hne connecting the target to the 
decay vertex ( P| ) of the four particles in the final state: 492 MeY/c < m < 504 MeV/c^ and < 100 {MeV/cf. 
After all the cuts we observed four events in the signal region. Figure || shows the m vs. the for events with 
all but these cuts. A linear extrapolation of the high data in this figure yields a background estimate of 0.25 ± 
0.10 events, consistent with the expectation from Monte Carlo studies. To further test the background estimate with 
higher statistics we removed the cluster shape cut and verified that the data matched the prediction in m and 
side bands. The probability of observing four events in the signal region due to fluctuation of the background is 
2.1 X 10~^, corresponding to a 4.2 a fluctuation of the estimated background. The branching ratio for Kl — > 
was calculated by normalizing the four signal events to a sample of Kl tt+tt^vt'^ events, collected with the prescaled 
normalization trigger. For the normalization events m-y^ was required to be within 3 MeV/c^ of m^ro, and the R||"^ 
and muon counter hit requirements were not enforced. The acceptance of these events was calculated to be 8.1% 
via Monte Carlo. We determined that (2.68 ± 0.04) x 10^^ Kl within an energy range of 20 to 220 GeV decayed 
between 90 and 160 meters from the target. The acceptance of the signal was (0.14 ± 0.01)%, so B{Kl fi'^fi^jj) 
= (10.4t^;^ (stat)) x 10"^ with m-y^ > 1 MeV/c^ which was the cutoff we used in generating the Monte Carlo events. 

We have calculated the branching ratio for this Kl Dalitz decay by performing a numerical integration of the 
tree-level {0{a)) Kl — > l^^^"^"^ matrix element with an ra^^ > 1 MeV/c^ cutoff. We performed a similar inte- 
gration of the Kl A*M7 matrix element, which included contributions due to virtual photon loops and emis- 
sion of soft bremsstahlung photons. Both integrations assumed unit form factors. The ratio of partial widths is 
2.789%. Multiplying this ratio with the measured value for B{Kl = (3-26 ± 0.28) x lO"'^ H yields 

B{Kl -> = (9.1 ± 0.8) X 10-^. 

The four-body phase space for Kl ji^ ^.^^^ can be parametrized by five variables, as in reference 

Figure ^ shows the distribution of the energy asymmetry of the photon pair (y-^), the angle between the normals 
to the planes containing the /i^/x^ and 77 in the center of mass (</>), and the minimum angle from any muon to any 
photon (Gmin)- The distribution of these kinematic variables for the four signal events is consistent with expectations. 

We examined several possible sources of systematic uncertainty in the measurement. The largest effects were due 
to a possible miscalibration of the calorimeter resulting in a mismeasurement of the photon energies, and particle 
identification. If we conservatively assume a 0.7% miscalibration of the calorimeter we obtain a 5.11% systematic 
error. The uncertainty due to muon identification was determined to be 4.2% by comparing the Kl flux with that 
obtained by using K^^ decays. The uncertainty in the Kl —> Tr+Tr^Tr" branching ratio is 1.59%. Adding these and 
other smaller contributions detailed in Table |l| in quadrature we assigned a total systematic uncertainty of 6.95% to 
the branching ratio measurement. 

In summary we have determined the branching ratio to be B{Kl n'^fi~"f"f) = (10.4t5;9 (stat) ±0.7 (sys)) x 10~^ 
with m^j > 1 MeV/c^. Defining the acceptance with a 10 MeV infrared cutoff for photon energies in the kaon 
frame {E*), our result is B{Kl ti+n--f-i]E* > 10 MeV) = (1.42tJ:° (stat) ± 0.10 (sys)) x 10"^. This is the first 
observation of this decay and is consistent with theoretical predictions. 
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TABLE I. The various backgrounds to Kl ji'^ fi, 77. n can can be mistaken for /u 
(P). Accidental clusters in the calorimeter identified as photons are designated ^acc- 



due to decay (D) or punch-through 



Decay 



Cause of fj, misid 



Events expected 



M^/U 77acc 

7r"'"7r~7r'' 



Kl 

f + — 
Kl ^ TT TT vr 

Kl — > n^iT~'K° 

Kl ^ ■k'^IJ^i' + 27acc 

Kl ■K^^'^V + 27acc 

Kl TT^TrV^i^ 
Kl tt^ttV^i^ 



DD 
DP 
PP 

D 

P 

D 

P 



0.093 ± 0.036 

< 0.056 

< 0.011 

< 0.011 
0.030±0.030 
0.032±0.032 

<0.005 
<0.004 



Total 



0.155 ± 0.081 



TABLE II. Systematic and statistical sources of uncertainty. Sources marked with (*) contribute to uncertainty in both the 
Kl flux and the acceptance for Kl — »■ fifi'yy relative to the acceptance for Kl — > tv^-k'tt^; other sources contribute only to the 
acceptance ratio. 



Source 
ZD 



Relative Uncertainty 



BKl — » tt" tt tt" 
Data statistics for Kl — > 7r"'"7r~7r° 
Simulation statistics 
Calorimeter scale and resolution 
Spectrometer scale and resolution 
Muon identification 
Signal trigger requirements 
Vertex quality requirement 
Spectrometer wire inefliciency 



1.59% (*) 
0.16% (*) 
0.22% (*) 
5.11% 
0.98% 
4.20% 
0.80% 
0.24% 
0.37% 



Total 



6.95% 
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FIG. 1. Energy deposited in the calorimeter by photons from Monte Carlo simulations of Kl 
accidental clusters (dashed) from data taken with a random trigger. 



fj.'^fj, 7 events (solid) vs. 
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FIG. 2. The distribution of Rjj^'^ (see text) from Monte Carlo simulations of the signal (sohd) and backgrounds from 
Kl — > 7r"'"7r~7r'' (dashed) that remain after the m~fj < 130MeV/c^ requirement. The arrow indicates the cut at -0.06, above 
which events were discarded. 




FIG. 3. The distribution of m^^ from Monte Carlo simulations of the signal (top) and backgrounds from /sT^s (bottom). The 
arrows indicate the cut at 340 MeV/c^, above which events were discarded. 
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FIG. 4. The cosine of the angle between the two photons, from Monte Carlo simulations of the signal (solid) and if,^3 
(dashed). The requirement of < 0.3 is indicated by the arrow. 
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FIG. 5. m vs. Px for the events that passed all other cuts. The box is drawn around the signal region, which contains 4 
events. 
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FIG. 6. The distributions for data (dots) and signal Monte Carlo (shaded) for the variables y^, Omin, and as defined in 
the text. 
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